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High resolution low-temperature absorption spectra of 0.2 % Pr'^"'' doped CsCdBrs were 
I measured in the spectral region 2000-7000 cm~^. Positions and widths of the crystal field 

levels within the '^-ffs, '^Ha, '^F2, and ^Fs multiplets of the Pr^+ main center have been 
determined. Hyperfine structure of several spectral lines has been found. Crystal field 
^ \ calculations were carried out in the framework of the semiphenomenological exchange charge 

QQ . model (ECM). Parameters of the ECM were determined by fitting to the measured total 

splittings of the '^Ha and ''Hq multiplets and to the observed in this work hyperfine splittings 
fT^ , of the crystal field levels. One- and two-phonon relaxation rates were calculated using 

' the phonon Green's functions of the perfect (CsCdBrs) and locally perturbed (impurity 

' dimer centers in CsCdBraiPr"^"'") crystal lattice. Comparison with the measured linewidths 

^— , confirmed an essential redistribution of the phonon density of states in CsCdBrs crystals 

' doped with rare-earth ions. 

> ' 
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^ ■ I. INTRODUCTION 

o 

Crystals of rare-earth (RE) doped quasi-one-dimcnsional double bromides CsCdBra are widely studied, 
mainly, because of their property to incorporate RE'^"'" ions in pairs, even at low RE concentrations. This 
makes them a promising material for up-conversion lasers. The structure of CsCdBra belongs to the 
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^ space group and consists of linear chains of face-sharing [CdBrg]^~ octahedra along the c-axis. Positional 
symmetry for Cd^+ is D-^d- RE'^+ ions substitute for Cd^"*" forming centers with different mechanisms 
of charge compensation. The main center consists of two RE'^^ ions placed in the chain on each side of 
an adjacent cadmium vacancy, [RE'^+-Cd^+ vacancy- RE3+] | |]. Both RE3+ ions in such a center are 
equivalent, their positional symmetry lowers from D^d to C'^y. Spectra of many other centers with different 
configurations were reported in EPR and optical studies of CsCdBraiRE crystals. In particular, positional 
symmetry for RE'^"*' ions in minor axial centers which are supposed to appear at higher RE concentrations, 
namely, [RE3+-RE3+-Cd2+ vacancy] or [RE^+- Gd^+ vacancy-Cd2+-RE3+], is again Cg^ but two R&+ 
positions are not equivalent, while in [RE^^-Cs"*" vacancy] centers positional symmetry of RE'^+ ion is Cs- 
Though the spectra of CsCdBr3:Pr'^+ were intensively studied before crystal field energies of the 

lowest excited multiplets were not investigated in absorption. The positions of crystal field levels for the 
main center, as reported by different authors, differ by 2-14 cm~^. Typical spectral resolution used in the 
mentioned works was 0.8-2 cm^^, while the only high resolution study (it concerned the level) revealed 
spectral lines as narrow as 0.1 cm^^ HiHl- Peculiar doublet line shape was observed for the '^H4{Ti) — > 
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"'^1)2 (ra) transition and explained by the combined effect of unresolved hyperfine structure and non- axial 
crystal strains 

The very first analysis of the crystal field energies of Pr3+ in CsCdBra H showed that the crystal field 
affecting the symmetric dimer centers consists of a strong cubic and a weak trigonal components and that 
the last one is determined mainly by the quadrupolar (Bj) term. This conclusion was confirmed in the 
studies of the crystal field energies of Tm'^^ and Ho'^+ symmetric pair centers where data obtained from 
both optical and submillimeter EPR investigations were taken into account . The sets of crystal field 
parameters obtained in Refs. |^,|| from the fitting of the simulated energy level schemes to the measured 40 
energy levels of Pr"^+ are close to one another and provide small enough root-mean square deviations (22.4- 
11.1 cm~^) between calculated and observed energy level data. However, some of the crystal field parameters 
in CsCdBr3:RE'^+ presented in the literature |^-^, in particular Bg and Bg, vary nonmonotonically with the 
occupation number of the 4/ electron shell. 

We have undertaken the low temperature high resolution infrared absorption study of CsCdBr3:Pr'^+ with 
the aims (i) to determine directly positions of crystal field levels in the region between 2000 and 7000 cm~^ 
and (ii) to study the line widths and shapes and thus to obtain information on relaxation processes and 
hyperfine interactions. We also performed crystal field calculations for the main center in the framework 
of the semiphenomenological exchange charge model (ECM) ||l0|, based on the analysis of the local lattice 
structure. Parameters of the ECM were corrected by fitting the calculated hyperfine splittings of crystal 
field levels to the experimental data from our measurements. The electron-phonon interaction effects, in 
particular, the low temperature relaxation of excited states of the (Pr'^+)2-dimer due to one- and two- 
phonon transitions between the crystal field sublevels of the '^i?5, ^Hq, and '^^3 multiplets, are studied 
in the framework of the recently derived rigid ion model of lattice dynamics of doped CsCdBraiRE crystals 




II. EXPERIMENT 

CsCdBr3 crystals containing 0.2% of Pr'^+ were grown by the method described in Rcf. Crystals 
easily cleave along the c -axis. We have prepared a 1.3 mm thick sample (sample A) with the c-axis parallel 
to the cleaved face. The sample was upheld at a controlled temperature of 5 K in a helium vapor cryostat. 
Unpolarized absorption spectra were measured with a BOMEM DAS. 002 Fourier transform spectrometer 
in the spectral range between 2000 and 7000 cm~^ with a resolution from 1.0 to 0.05 cm~^. In Fourier 
transform spectroscopy, the quantity Res= 1//3, where C is the maximum optical path difference, is indi- 
cated as the resolution. The full width at half hight (FWHH) of the instrumental function depends on a 
particular apodization function used. We used no apodization and, thus, worked with the narrowest possible 
instrumental function (FWHH= 0.6 Res). It is easy to show that, in this case, the shape of a spectral line 
is essentially unchanged provided that FWHH < Res. 

Precision of the experimental line positions was 0.05-0.3 cm~^, depending on a particular line. Precise 
absolute wave number scale is an intrinsic property of Fourier transform spectroscopy. Line widths were 
determined from absorbance spectra calculated with a zero-line taken at a half of transmitted intensity. In 
such a way we tried to take into consideration the fact that the incident light is unpolarized but a majority 
of lines is 100% polarized. 

The lowest frequency region of the ^H^ — > ^F^ spectral transition has also been measured with the 
resolution of 0.005 cm^^ using another sample, 3.1 mm thick and oriented approximately along the c-axis 
(sample B). It allowed almost pure k||c geometry and, consequently, E, H _L c polarization of the incident 
light. 



III. EXPERIMENTAL RESULTS 

Fig. 1 shows low temperature transmittance spectra of CsCdBr3:Pr'^+ (0.2%) corresponding to different 
infrared transitions. The first excited state of the ground multiplet lies at about 170 cm~^ a-nd is 
not populated at helium temperatures. So, low temperature absorption spectra display crystal field levels of 
the excited multiplets directly. They are listed in Table 1, together with data from other publications. 
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As the ground state is the ^H4{^Ti) singlet, excited Fi states manifest themselves in 7r-polarization but 
Fa states — in cr-polarization, while F2 states are silent in electric dipole approximation (see Table 2). 
Magnetic dipole transitions are usually very weak, with a possible exception for the '^H4 — > ^H^ transition 
that is allowed for the free Pr'^+ ion. 

The symmetries of the crystal field levels have been determined from a comparison between the spectra 
of the two differently oriented samples (samples A and B). The assignement of the 6483.3 cm~^ level as a F3 
level is strongly supported by the well resolved hyperfine structure (hfs) observed in high resolution spectra 
(see inset 1 of Fig. Id and Fig. 2a). The point is that while F3 levels exhibit the magnetic hfs, Fi and F2 
do not. The next F3 level also shows the traces of unresolved hfs, namely, the 6503.5 cm"^ line has a flat 
bottom (see inset 2 of Fig. Id). 

Expansion of the 6483.3 cm^^ line into six components yields 0.040 cm^^ wide Gaussians. The Gaussian 
shape of hfs components gives evidence that an inhomogeneous broadening due to random crystal fields 
exceeds a homogeneous broadening due to nonradiative transitions to the nearest ^F2 level lying at about 
10 cm~^ lower ||,|^ (we did not observe '^i?4(^Fi) '^i^3(F2) transitions forbidden as electric dipole ones). 
The lineshapes of the 6496.5 cm~^ (Fi — > Fi) and 6503.5 cm~^ (Fi — > F3) lines reveal homogeneous broaden- 
ing of 0.17 and 0.31 cm^^ , respectively. Such a broadening comes, evidently, from electron-phonon interaction 
and implies a noticeable density of phonon states in the region between 20 and 30 cm~^ corresponding to 
the distances from the considered levels to lower crystal field levels. 

Other examples of level broadening due to electron-phonon interaction can be found in the '^i?5 multiplet 
which is situated in the region of about 2400 cm~^ (see Fig. la). The lowest level of this multiplet lying 
at about 2234 cm~^ was reported to be the ^F2 level, from polarized fluorescence measurements [0. In 
that case, a strong very narrow (with the width of 0.05 cm^^) line observed at 2234.8 cm~^ corresponds 
to the magnetic dipole transition ^H^i^Ti) ^Hr,{^r2). As we have already mentioned, such a situation 
is possible just for the '^Hi '^H^ spectral multiplet. The next level lies at 26 cm~^ from the bottom 
of the multiplet and, though suffers from a phonon relaxation, exhibits hfs as it follows from the shape of 
2261.4 cm~^ line. This is in accordance with the previously made assignments of this level symmetry as 
F3. Three closely spaced levels at 2316.6, 2331.9 and 2347.8 cm^^ are broad. Their widths are 3.0, 3.3 and 
3.7 cm^^ respectively, and it is not possible to determine their symmetries. Calculations described below 
show that an overall hyperfine splittings of F3 levels do not exceed 0.5 cm~^. Consequently, the level widths 
of 3-4 cm~^ are, mainly, due to rapid phonon-assisted decay to lower levels. In accordance with the crystal 
field calculations (see below) , at least one of these three broad absorption lines is to be assigned to another 
center. The next level is more narrow again. According to the shape of the line at 2546.7 cm~^ it contains 
unresolved hfs and thus corresponds to the F3 level. The distances from this level to lower levels exceed the 
length of the phonon spectrum (~ 180 cm~^ [01)- Because of that, one-phonon relaxation is not possible. 

The highest frequency line observed at 2572.1 cm~^ could be attributed to the '^i?4(^Fi) — > '^H^{^T2) 
magnetic dipole transition. But the ^iJ5(^F2) level lying at 26 cm^^ above the '^ff5(^F3) level (2546 cm^^) 
should be broadened due to one- and two-phonon transitions to the lower crystal field sublevels, and the 
width of the respective spectral line should be about 0.3-0.5 cm^^ (see Table 3), not 0.05 cm~^ as observed. 
Probably, this line belongs to some other Pr'^+ center. Weak narrow lines observed at 2153.6 and 2167.3 cm~^ 
(see Fig. la) also have to be assigned to some other centers. If these levels belonged to the main center, the 
level at 2234.8 cm~^ would be broadened by the nonradiative transitions to these levels. It is not the case 
experimentally. 

We were able to find only two levels in the multiplet and the lowest sublevel in each of "^^4 and '^Hq. 
The level at 5073 cm~^ which lies far (250 cm~^) from a lower level is, evidently, broadened by unresolved 
hfs and is thus a F3 level. The measured hyperflne splittings and widths of hyperfine sublevels found from 
the experimental line shapes are presented in Table 3, to compare with the calculated hyperfine splittings 
and the estimated one-phonon decay rates (see the next section). 

Many of the observed levels coincide with those of Refs. found from selectively excited fluorescence 
measurements. We did not observe any additional structure due to noncquivalence of Pr'^+ positions in 
asymmetric pair centers. Thus, our results support the conclusion of Refs. jfj-^ that the main center in 
CsCdBr3:Pr'^+ is a symmetric pair [Pr'^+~Cd^+ vacancy-Pr'^+] . 
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IV. CRYSTAL FIELD CALCULATIONS 



The energy level pattern of the i4ipj.3+ j-^j^g nuclear spin / — 5/2, 100% abundant, the 4/^ electron 
shell) in the trigonal crystal field can be represented by eigenvalues of the effective parametrized Hamiltonian 



H = Vee+ Hso + ^E{L, S, J) + H,f + Hhf, 



(1) 



where T4e is the two-body Coulomb energy; Hso corresponds to the spin-orbit interaction (we use in the 
present work ( = 746.2 cm^^ and Slater integrals F2 = 304.4; F4 = 45.47; Fq = 4.41 cm~^ Q); H^f is the 
crystal field Hamiltonian 
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(Op are the Stevens operators ||l0[| , the quantization axis z coincides with the symmetry axis C3 of the crystal 
lattice), and Hhf is the nuclear energy responsible for the hyperfine structure of the optical spectra. Because 
we neglect many minor interactions (three-body terms, spin-other orbit interaction, two-body correlations 
in the crystal field terms and etc.), shifts of the multiplet centers of gravity AE{L,S,J ) are introduced 
to fit the calculated spectrum to the experimental data. Only projections of the electron-nuclear magnetic 
and quadrupolar interactions on multiplet manifolds with fixed orbital (L), spin (S) and total (J) angular 
moments were considered: 



Hhf = A{L,S, J)JI 



e^Q(l~7) 
4/(2/- 1) 



i^K, [3/2 -/(/+1) 



4/(2/ - 1) 



I (3 Jf - J( J + 1)) (342 - /(/ + 1)) + 1 {JlH + Jill) 

i (J, j+ + j+j,) (/,/_ + /_/,) + i (J, j_ + j_ J,) (/,/+ + /+/,; 



(3) 



Here A{L,S,J) is the magnetic dipole hyperfine constant (in the free ion A{^H4) — 1.093 GHz [|lj], the 
corresponding hyperfine constants for other multiplets are obtained using this value and neglectin g th e 
core electron polarization contributions which do not exceed 1.5% of the 4/ electron contributions fl^]), 
Q = —5.9 X 10""^" m^ |l^ is the nuclear quadrupole moment, 7 = — 70± 10 is the Sternheimer antishielding 
factor |l|], aj are reduced matrix elements of the second rank spherical operators, (r^) are moments of the 
radial wave function of the 4/ electrons calculated in |17|, and the electric field gradient at the Pr^+ nucleus 
equals 



3cos^e{i) - 1 
S(ij3 ^ 



(4) 



where q{i) is the charge (in units of the proton charge e) of a lattice ion with spherical coordinates R{i), 0{i) 
and (p{i) in the system of coordinates having its origin at the Pr'^+ nucleus. 

As we have already mentioned in the Introduction, the crystal field parameters obtained earlier from the 
fitting to the measured energy levels of Pr'^+ |^,|| and Nd^+ |^ in CsCdBrs, though providing small enough 
deviations between the calculated and observed energy levels, exhibit nonmonotonic variation within the 
RE'^+ family (see Table 4). Moreover, there are large discrepancies between the measured in this work hfs of 
the absorption lines and the hfs, calculated with the crystal field parameters from Ref. In particular, the 
6483 cm~^ line exhibits well resolved hfs of total width 0.21 cm~^, while the calculated width is 0.11 cm~^. 
The formal fitting procedure employed in |^,^j9j evidently brings about overestimated 6-rank parameters of 
the trigonal crystal field component. 

To clear up relations between the crystal field parameters and their origin, we estimated them in the 
framework of the exchange charge model |l^. The crystal field is represented as a sum {Bp = B^^ 
of the electrostatic field of the lattice ions and the exchange charge field defined by the parameters 
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and 

2(2p+l) 2;^fc^ 2 , ^ ..^2 , ^ o ,.^2^ Qp(^W'V^(^)) 

5pS = e Kp2_^{G,Ss[i) +GaSa{i) +jpG^S^{i) ) — , (6) 

i 

respectively. Here are numerical factors [|lO|, ap are the shielding factors, 72 = 3/2, 74 = 1/3, 76 = —3/2, 
and Go- = Gjr = Gp are dimensionless parameters of the model. The expression (6) involves the sum over 
only the nearest neighbors (six bromine ions) of the RE^+ ion. The overlap integrals Ss = {4:f,iTi = 0|4s), 
Sa — {4:f,m = 0|4p, m — 0), and Syr — (4/, m — l\Ap,m = 1) have been computed using the radial 4/ — 
wave function of the Pr'^+ ion from Ref. |17| and the 4s, Ap wave functions of the Br~ ion given in p^ . 



Dependencies of the overlap integrals on the intcrionic distance R (in Angstroms) are approximated by the 
following expressions: 

Ss = 0.09988 exp(-0.20882i?2i7238). 

S„ = 0.04610 exp(-0.03733i?3°662). (7) 
= 0.39635 exp(-0.87514i?i-30253)^ 

Equations (H, |^, |^) present the crystal field parameters as explicit functions of relative positions of the 
impurity Pr'^+ ion and the lattice ions. We analyzed distortions in the relaxed CsCdBrs lattice containing 
the substitutional dimerfPr'^^-Cd^^ vacancy-Pr'^+] in the framework of the same quasimolccular model as 
the one already used in |^| to calculate the structure of symmetric Tm^+ dimers. Values of ion displacements 
from their equilibrium positions in the perfect lattice were obtained by minimizing the potential energy of 
the cluster which involved the nearest 50 ions around the symmetric dimer center. These ions belong to the 
first (12 Br^), second (2 Cd^^), third (12 Cs+), and fourth (24 Br~) coordination shells of the impurity 
Pr'^+ ions (see Table 5). 

The energy of the intcrionic interaction is assumed to be a sum of Coulomb and non-Coulomb terms. The 
effective ion charges {q{C(l^+)= 1.56, q (Cs+)= 0.84, (?(Br-)= -0.8, q{Vr^+)^ 2.34) and the non-Coulomb 
force constants, corresponding to interactions between the nearest neighbors, Pr^+-Br~, Cd^^-Br", Cs"*"- 
Br~, and Br~-Br~, have been obtained from the study of the CsCdBrs lattice dynamics |jll|. According 
to results of simulating the relaxation energy minimum with respect to the ion displacements (see Table 5), 
the distance between the Pr^+ ions attracted by the Cd^^ vacancy diminishes from the value of the lattice 
constant c = 0.6722 nm ||l^ down to 0.5915 nm. This value is very close to the distance of 0.593 nm between 
the Gd'^"'' ions and 0.5943 nm between the Tm'^+ ions 1^,^ in the symmetric pair centers in CsCdBra, 
as determined from the EPR spectra. In the perfect lattice, the first coordination shell of a Cd^"*" ion has 
the radius of 0.277 nm; in the symmetric dimer center six Br~ ions nearest to the Cd^^ vacancy are forced 
outwards, but the displacements of the Pr^+ ions along the dimer axis are so large that the corresponding 
Pr'^+-Br~ distance decreases down to 0.275 nm. On the contrary, the terminal triangles of the Br~ ions lag 
behind the impurity ions and the corresponding intcrionic distance increases up to 0.290 nm. These values 
agree with the distances (0.264 and 0.286 nm) between the Yb^"*" and the nearest Br^ ions in the intrinsic 
dimer units in Cs3Yb2Br9 

Taking into account the local lattice deformation, we obtained from Eqns. (|[ |[ 0) the following crystal- 
field parameters (values are in cm~^) for Pr ion in the symmetric dimer center (the electrostatic contri- 
bution to the quadrupolar crystal field component = —6^(1 — (J2){t'^)VzzIA was computed exactly with 
the Ewald method): 

B° = -445.5(1 - (72) - 1.123Gs - 2.828Gp; 
Bl = -42.7(1 - 0-4) - 4.52G^ - 10.12Gp; 
Bl = 1184(1 - (74) + 156.84G^ + 341.60Gp; 
Bl = 1.68(1 - CTe) + IMGs + 0.4395Gp; 
Bl = 13.4(1 - o-g) + 12.56G, + 13.56Gp; 
Bl = 24.9(1 - <7g) + 27.61Gs + 9.06Gp. 

Values of the model parameters Gs — 5.4; Gp = 6.1 have been determined by comparing the calculated and 
measured total splittings of the ground state '^7?4 and of the excited multiplet '^Hq, the shielding factors 
have been fixed at (T4 = (Tg = 0, (T2 = 0.89 in accordance with the theoretical estimations |16|. Thus 



5 



obtained crystal field parameters (column a of Table 4) were used as starting values for the fitting procedure 
constrained by the additional conditions that the calculated hfs had to agree with the experimental data. 

The final values of the crystal field parameters for the Pr^+ centers as compared with the corresponding 
parameters of Tm'^+j Ho^+ and Nd'^"'" centers are given in column h of Table 4. Despite the large trigonal 
distortion of the nearest bromine octahedron, the main features of the crystal field splittings are determined 
by the cubic component of the crystal field. The Sg and parameters involve the largest changes from 
the initial (calculated) to the final (fitted) values. These differences are most probably caused by the 
overestimated displacements of the six Br7-Bri2 ions in the basis plane (see Table 5). The calculated energy 
levels of the Pr'^+ ions presented in Table 1 were obtained by diagonalizing the Hamiltonian (1) in the space 
of 516 states (terms , ^F, ^P, ^G, were taken into account) with fitted barycenters of the free-ion 
multiplets. The measured splittings of the ^H^, ^H^, ^Hq, ^Fj, ^^3, ^^4, ^Pi, ^Pj, multiplcts of Pr3+ 
are satisfactorily described by our final set of crystal-field parameters (column b in Table 4). 

The computed hfs of Fi and r2 crystal field singlets consists of three doublets with the total width AE of 
no more than 0.02 cm~^ which can not be resolved in the optical spectra. The quadrupolar contributions to 
the hfs are very small (~ 10"'* cm"-'^ ), the calculated hfs of Fa doublets is actually equidistant with intervals 
of no more than 0.1 cm~^. An example of the simulated envelope of the absorption line ■^i?4(^Fi) ■^F3(^F3), 
where individual transitions between different electron-nuclear states are presented by Gaussians of 0.04 cm~^ 
width, is given in Fig. 2b. The observed irregular shape of the hfs of some absorption lines (see Fig. 2a and 
Ref . [gl ) is most probably caused by crystal fields of low symmetry. We discussed this point in more details 
in another publication | p3[ |. 

The most essential difference between results of our calculations and fitting procedures performed in |^,|| is 
connected with assignments of symmetry types of some crystal field energy levels. In particular, symmetries 
of the third {^T^} and fourth (Fi) sublevels of the ^H^ multiplet and of the second (^F3) and third (F2) 
sublevels of the ground "^H^ multiplet are inverted in |^,|| (it should be noted that our scheme of crystal field 
energies of the ground state agrees with the assignment of Ref. ) . Additional measurements of polarized 
spectra are to be carried out to check our predictions. 



V. RELAXATION BROADENING OF THE CRYSTAL FIELD LEVELS 



Even at the liquid helium temperature the hfs of the most of excited crystal field sublevels is masked by the 
spontaneous relaxation broadening. Strong electron-phonon interaction effects in CsCdBr3:RE'^+ crystals 
originate from the specific density of phonon states which has large maxima in the low frequency region (20- 
40 cm^^) in the perfect crystal lattice (see Fig. 3 and Ref. ||ll[l). The Hamiltonian of the electron-phonon 
interaction, expanded in a power series in ion displacements from their equilibrium positions, can be written 
as follows 

-ffel-ph = Y^Va{s)[Ucis) - Ua(RE)] + 

sa 

+ ^E^"/5(^) - "a(RE)][u^(s) - U0(RE)]...; (8) 

sq/3 

pk pk 

where [u(s) — u(RE)] is the difference between dynamic displacements of the ligand ion s and the RE ion, 
and Bp^^{s), ^^{s) are the coupling constants. The probability of the one-phonon transition between the 
initial (i) and final (/) states of the RE ion with the energy gap hujif > can be represented as 

= f E {f\Vo.{s)mm9Mss'\u;.f){t\Vf,is')\f){l + n{u;,f)), (9) 

sas' P 

where n(a;) is the phonon occupation number, and gap(ss'\Lu) are spectral representations of the Green's 
functions for differences between the ion displacements |24|. We performed a calculation of the relaxation 
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rates for all the crystal field sublevels within the manifolds of '^i?5, ^-^2 and ^F^ multiplets using the phonon 
Green's functions of the perfect (CsCdBra) and locally perturbed (impurity dimer centers in CsCdBr3:Pr'^+) 
crystal lattices obtained in . The formation of a dimer leads to a strong perturbation of the crystal lattice 
(mass defects in the three adjacent Cd^+ sites and large changes of force constants). As it has been shown 
in , the local spectral density of phonon states essentially redistributes and several localized modes corre- 
sponding to different representations of the dimer point symmetry group D^d appear near the boundary of the 
continuous phonon spectrum of the unperturbed lattice. As an example of changes in the phonon spectrum, 
we present in Fig. 3 spectral densities Daai'-^) of displacement-displacement autocorrelation functions 

{uliR)) = J Dc4cj)duj, (10) 

in the perfect {R — Cd^"*") and perturbed {R = Pr^+) lattices at zero temperature. The contributions to 
Dxx{i^) from the localized modes of T^g (196.5 cm~^) and r3„(190 cm~^) symmetry, and to Dzz{uj) from 
the localized modes of Fig (191.7 cm~^) and r2ti(191.2 cm~^) symmetry are presented by Lorentzians with 
the proper weights and the full width of 2 cm^^. 

We took into account interactions of Pr'^+ ion with its nearest neighbors (six Br~ ions) only. Values 
of coupling constants were computed in the framework of the ECM with the same parameters as those 
used in the crystal field calculations (explicit expressions of the Bp^{s), ^^^{s) are given in @j24)). 
Matrix elements of electron operators Va{s) were calculated with the eigenfunctions of the Hamiltonian 
(1). Calculated inverse lifetimes l/n — Wi = J^f^if (the sum is over states / which belong to the 
same multiplet as the crystal field level i does, with lower energies) are given in Table 3. It is seen that 
despite many simplifying approximations there is a good correlation between the measured linewidths and 
the relaxation rates obtained with the perturbed Green's functions. Calculations with Green's functions of 
the perfect lattice gave overestimated (up to an order of magnitude higher) values of relaxation rates. Thus, 
we conclude that the increased Pr''+-ligand elastic interaction, as compared to the Cd^+-Br~ interaction, 
and corresponding enhancement of correlations between displacements of the impurity RE'^+ ion and its 
neighbors strongly suppresses the electron-phonon coupling. 

It should be noted that the measured width of the doublet ^r3{^H^) is five times larger than the estimated 
total hfs width, though the one-phonon relaxation broadening of this level is not possible at low temperatures. 
We suppose that this level, and two upper levels ^r2('^7?5), '^T3{^Hr,) are essentially broadened due to the 
two-phonon relaxation. The probability of the spontaneous emission of two phonons contains three terms 
which correspond to the first order contribution from the nonlinear electron-phonon interaction: 
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lmga-f{ss' \uji)lmgi3s{ss' \uj2)S{uJif - loi - W2)dwidw2, (11) 
to the second order contribution from the linear electron-phonon interaction: 

J J \ {oJij - UJl){uJu - Ul) {cJij - UJl){u}a - CJ2)y 

d{ujif — uji — uj2)dujiduj2, (12) 

where 

Mf^M = E {f\Va{s)\j)lmg^0{ss'\uj^){l\V0{s')\m), (13) 
and to the combined action of the first and second order transition amplitudes: 

= ^ E E((/i^"M«)N)(*in-(s")U)(j|v^7(^')i/)+c.c.)x 

I [ lmgaj{ss'\uJi)lmgi3s{ss"\Lj2)-. -S{uJif - uJi - U!2)di^iduj2- (14) 



7 



The calculated sum of two-phonon relaxation rates X]/ ^i/ corresponding to spontaneous transitions 
from the doublet ^T^^'^H^) to all lower sublevels of the '^H^ multiplet, induced by the nonlinear terms in 
Eqn. (H) equals 1.8 x 10^" s^^ and is almost exactly cancelled by the contribution —1.7 x 10^° s^^ from 
" crossed" terms of Eqn. (|4|) . We have obtained the remarkably larger contribution to the inverse lifetime of 
this doublet (7.1 x 10^" s ), that agrees satisfactorily with its measured linewidth, from the second order 
terms given in Eqn. (p^). It should be noted that relaxation processes that involve excitations of localized 
modes of and T^g symmetry dominate yielding more than a half of the total relaxation rate. 

VI. CONCLUSION 

High resolution (1-0.005 cm^^) infrared (2000-7000 cm-^) absorption spectra of CsCdBr3:Pr3+ (0.2%) 
were taken at low temperatures. Crystal field levels of the '^-ffs, ^Hq, ^F2, and ^F^ multiplets for the 
symmetric pair center [Pr'^"'"-Cd^"'" vacancy-Pr'^+] were found directly in absorption. Hyperfine structure, 
inhomogeneous and relaxation-induced widths of several crystal field levels were measured. 

We performed crystal field calculations in the framework of the semiphenomenological exchange charge 
model and took into account local lattice deformation around an impurity center. Only two parameters of the 
model had to be determined by a comparison of calculated and experimental energy levels. Thus obtained 
crystal field parameters were corrected by fitting the calculated hyperfine structure to the measured one. 
It follows from our results that despite a large trigonal distortion of the bromine octahedra in the nearest 
surroundings of impurity Pr^+ ions, crystal field splittings are determined, mainly, by the cubic component 
of the crystal field. 

Using the results of the crystal field calculations and of the lattice dynamics analysis we were able, for the 
first time, to calculate exactly (without any additional fitting parameters) one- and two-phonon spontaneous 
relaxation rates. The calculated total relaxation rates reveal the linewidths which are in reasonable agreement 
with the measured ones. The obtained information on relative efficiency of linear and nonlinear terms in the 
Hamiltonian of the electron-phonon interaction in stimulation of the phonon emission is important for the 
thorough derivation of the theory of multiphonon relaxation in rare-earth compounds. 

The peculiarities of the electron-phonon interaction effects in CsCdBr3:Pr^+ crystals originate from the 
specific density of phonon states that extends by only ~ 180 cm~^ in the perfect CsCdBra lattice, and from 
the localized modes induced by the impurity dimer centers. Another important feature and, probably, the 
most interesting result of this study, is a strong suppression of the effective electron-phonon coupling due to 
a local increase of elastic forces in the activated crystal and the corresponding enhancement of correlations 
between displacements of the impurity RE ion and its neighbors. 

ACKNOWLEDGEMENTS 

An encouraging interest of G. N. Zhizhin is acknowledged. This work was supported in part by the RFBR 
(Grant N 99-02 16881) and the Russian Ministry of Science (Grant N 08.02.28 of the Program "Fundamental 
Spectroscopy"). 



[1] F. Ramaz, R. M. Macfarlane, J. C. Vial, J. P. Chaminade, and F. Madeore, J. Lumin. 55, 173 (1993). 

[2] E. Antic-Fidancev, M. Lemaitre-Blase, J. P. Chaminade, and P. Porcher, J. Alloys Compd. 225, 95 (1995). 

[3] K. M. Murdoch and N. J. Cockroft, Phys. Rev. B 54, 4589 (1996). 

[4] J. Neukum, N. Bodenschatz, and J. Heber, Phys. Rev. B 50, 3536 (1994). 

[5] J. P. Chaminade, R. M. Macfarlane, F. Ramaz, and J. C. Vial, J. Lumin. 48 & 49, 531 (1991). 

[6] J. Heber, M. Lange, M. Altwein, B. Z. Malkin, and M. P. Rodionova, J. Alloys Compd. 275-277, 181 (1998). 



8 



[7] B. Z. Malkin, A. I. Iskhakova, V. F. Tarasov, G. S. Shakurov, J. Heber, and M. Altwein, J. Alloys Compd. 

275-277, 209 (1998). 
[8] G. W. Burdick and F. S. Richardson, J. Alloys Compd. 275-277, 379 (1998). 

[9] J. R. Quagliano, N. J. Cockroft, K. E. Gunde, and F. S. Richardson, J. Chem. Phys. 105, 9812 (1996). 
[10] B. Z. Malkin, in Spectroscopy of Solids Containing Rare-Earth Ions, edited by A. A. Kaplyanskii and R. M. 

Macfarlane (North-Holland, Amsterdam, 1987), p. 13. 
[11] B. Z. Malkin, A. I. Iskhakova, S. Kamba, J. Heber, M. Altwein, G. Schaack, submitted to Phys. Rev. B. 
[12] R. B. Barthem, R. Buisson, F. Madeores, J. C. Vial, and J. P. Chaminade, J. Phys. (Paris) 48, 379 (1987). 
[13] G. Gorller-Walrand and K. Binnemans, in Handbook on the Physics and Chemistry of the Rare-Earths, edited 

by K. A. Gschneidner and LeRoy Eyring (North-Holland, Amsterdam, 1996), Vol. 23, p. 121. 
[14] M. A. H. McCausland and I. S. Mackenzie, Adv. Phys. 28, 305 (1979). 
[15] G. Netz, Z. Phys. B 63, 343 (1986). 

[16] V. C. Das and A. K. Ray Chaudhuri, J. Phys. C 6, 1385 (1973). 
[17] A. J. Freeman and R. E. Watson, Phys. Rev. 127, 2058 (1962). 

[18] E. Clementi and C. Roetti, Atomic Data and Nuclear Data Tables (Academic Press, New York, 1974), Vol. 14, 

N 3-4. 

[19] G. L. McPherson, A. M. McPherson, and J. A. Atwood, J. Phys. Chem. Sol. 41, 495 (1980). 
[20] L. M. Henhng and G. L. McPherson, Phys. Rev. B16, 4756 (1977). 

[21] V. F. Tarasov, G. S. Shakurov, B. Z. Malkin, A. I. Iskhakova, J. Heber, and M. Altwein, JETP Lett. 65, 559 

(1997). 

[22] M. P. Hehlen and H. U. Gudel, J. Chem. Phys. 98, 1768 (1993). 

[23] E. P. Chukalina, M. N. Popova, E. Antic-Fidancev, and J. P. Chaminade, Phys. Lett. A, 258, 375 (1999). 
[24] M. N. Popova, E. P. Chukahna, B. Z. Malkin, S. K. Saikin, Phys. Rev. B 61, 7421 (2000). 



9 



FIGURES 




2234 2235 



(a) 

I I I I I I I I I I I I I I I I I I I I I I I I I I I 

2200 2300 2400 2500 2600 
Wave number (cm"^) 



10 



U 0.80' 
H 



CZ3 

f5 



CZ3 



(b) 



4400 



4500 



4600 



0.80 



0.40- 



(c) 



5100 



5150 



I 0-80' 



0.00 



3tT 3t:' 
H4 ^ F3 



w , 




65036504 



6482.8 6483.6 



(d) 



T" 



6500 



6550 



Wave number (cm" ) 



FIG. 1. Transmittance spectra of CsCdBrsiPr^"'" at 5K corresponding to the optical transitions from the ground 

state ''H4(^ri) to the excited (a) ^Hs, (b) ^Hg, (c) ^F2, and (d) ''Fa crystal field multiplets. Spectral resolution is 
0.055, 0.3, 0.3, and 0.05 cm~^ for (a), (b), (c), and (d), respectively, and 0.005 cm^^ for insert 1 of Fig. 1(d). Very 
sharp lines in the spectral region 2300-2400 cm~^ are due to residual CO2 in the spectrometer. 
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FIG. 2. (a) Measured with resolution 0.005 cm ^ and (b) calculated hyperfine structure of the 
^HiCVi) F3CT3) transition of CsCdBr3:Pr3+. 
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FIG. 3. Simulated spectral densities of displacement-displacement autocorrelation functions for Cd^"*" in CsCdBrs 
(dotted curves) and for Pr^+ in the impurity dimer centers (solid curves). 
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TABLE I. Positions (cm ^) and symmetries of energy levels for the symmetric dimer centers, and widths (in 
brackets, cm~^) of absorption lines from the ground state in CsCdBrsiPr^"'". 









Measured 








Multiplct 


This work 




Literature data [1-4] 


Calculated 




'Hi 










Ti 





^Fi 








166-170 


r2, Ta ? 


168 


'Fa 








194-196 


Fa, r2 ? 


188 


F2 








323-328 


La 


293 


'Fa 








551-555 


Ti 


545 


'Fi 








574-580 


Fa 


564 


'Fa 




2234.8 (0.05) 




? 2230-2236 


F2 


2235 


'F2 




2261.4 (0.64) 




2256-2264 


Fa 


2251 


'Fa 




2316.5 (3.0) 




? 2309-2318 


Fi, Fa? 


2330 


'Fa 




2331.9 (3.3) 


Fa, Ti 


? 2327-2334 


Fa, Fi? 


2334 


Fi 




2348.4 (3.7) 




2347 










2546.7 (0.7) 


Ts 


2533-2547 


Fa 


2540 


'Fa 








2544-2595 


F2 


2593 
2620 


'F2 
"Fa 




4362.8 (0.19) 


Fa 


4362-4364 


Fa 


4363 


'Fa 








4374-4378 


Li 


4390 


'Fi 








4413-4415 


Fa 


4414 


'Fa 








4494 


F2 


4539 


'F2 








4700-4704 


Fa 


4710 


'Fa 








4725-4730 


Ti 


4733 


'Fi 








4761 


Fa 


4765 


*Fa 








4781-4808 


F2 


4804 


'F2 








4808-4843 


Li 


4815 


'Fi 
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Ta 
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Fa 


5073 


'Fa 
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Ti 
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'F2 
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Ta 
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Fa 
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Ti 
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Fi 
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6526-6535 


F2 
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'F2 
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Fi 
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'Fi 
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F2 


6907 


F2 








6903-6908 


Fa 


6902 


'Fa 








6914 


Fa 


6930 


'Fa 








7103 7104 


Ti 


7115 


'Fi 








7111-7117 


Fa 


7108 


'Fa 


'G4 






9777 




9777 
9869 
9899 

9955 

10427 

10453 


'Fi 
'Fa 
F2 

'Fa 
'Fi 
'Fa 








16536-16540 


Fa 


16536 


'Fa 








16567-16570 


Fi 


16588 


Fi 








17004-17011 


Fa 


17001 


'Fa 








20386-20393 


Fi 






'Pi 






20956-20964 


F2 


20964 


F2 








21011-21047 


Fa 


21010 


Fa 


'P2 






22106-22118 


Fa 


22117 


'Fa 








22163-22177 


Fi 


22165 


Fi 








22216-22229 


Fa 


22223 


2p 
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T'ARTiT^l TT Splpptinn mips for plprtTir Hinnlp 


(Hi anH ma.p'npt.ip Hiriolp (//I t.ransit.ioTiR in 
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TABLE III. The calculated and measured total hyperfine structure widths AE/?i, the one-phonon relaxation rates 
W and the widths Act) (10^" s^^) of the crystal field sublevels at the temperature of 5 K. 

^^+'Lj Crystal field AE//i W^"^ Aw 

energy (cm~^) Calculated Measured Calculated Measured 



^Ea (2235) 


0.05 


<0.94 





0.94^"^ 


^Fs 26 


6.68 


6.41 


15.5 [113] 


5.65 


^Fs 82 


5.43 




255 [573] 


56.6 


^Fi 97 


0.16 




157 [745] 


62.2 


^Fa 312 


2.65 




72W 


13.2 


^Fa 358 







52 [59] 




"Eg 385 


9.03 




59 [311] 




^Eg (4363) 


2.00 


2.26 





LSO^*"^ 


^Ea (5073) 


1.15 


1.32 





0.94(''^ 


^Eg 75 


2.09 




14 [31] 




^Ei 80 


0.03 




32 [57] 


52.8 


^Ea (6473) 


0.03 









^Eg 10 


3.71 


3.96 


0.52 [9.0] 


0.75(''^ 


^Ei 23 







4.5 [26] 


3.20 


^Eg 30 


1.53 


1.69 


15.8 [95] 


5.85 


^Ea 53 







63 [98] 





- numbers in square brackets were obtained with the phonon spectrum of the unperturbed crystal lattice, 

- the inhomogeneous width, 
('^^ - two-phonon relaxation rate. 



TABLE IV. Crystal field parameters Bp (cm '^) of symmetric dimer centers in CsCdBrg rRE^"*" crystals. 



p 


k 


Tm (4/12) 
[6] 


Ho (4/1°) 
[7] 


Nd (4/3) 
[9] 


[8] 


Fr (4/2) 

a 


this work 

b 


2 





-81.4 


-75.6 


-102 


-79 


-72.3 


-70 


4 





-82.0 


-94.3 


-121 


-149 


-128.S 


i -140.6 


4 


3 


2362 


2841 


3698 


4044 


4115 


4191 


6 





11.26 


12.75 


14.06 


22.62 


14.95 


15.1 


6 


3 


153 


216 


342 


347 


164 


279 


6 


6 


106 


150 


187 


37 


229 


101 
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TABLE V. Structure of the symmetric dimer Pr'^"^-Cd^^ vacancy-Pr^^ (lattice constants a — 0.7675 nm; 
c = 0.6722 nm; v = 0.1656 ||l|l). 



Ion 


X j Cb 


11 1 n 

y / ^ 


z 1 c 
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+ [05 — 060031 
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+a/3[1 /a/3 + 090511/9 
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(Cs+ 
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(Br" 
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